Properties of charmonia in a hot equilibrated medium 
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We investigate the properties of charmonia in a thermal medium, showing that with increasing 
temperature the decay widths of these mesons behave in a non-trivial way. Our analysis is based 
on a potential model with interaction potential extracted from thermal lattice QCD calculations of 
the free-energy of a static quark-antiquark pair. We find that in the crossover region some decay 
widths are extremely enhanced. In particular, at temperatures T ~ T c the decay widths of the J/9 
that depend on the value of the wave function at the origin are enhanced with respect to the values 
in vacuum by about a factor 2. In the same temperature range the decay width of the process 
Xcj — * J/9 + 7 is enhanced by approximately a factor 6 with respect to the value in vacuum. At 
higher temperatures the charmonia states dissociate and the widths of both decay processes become 
O |' vanishing small. 
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O I. INTRODUCTION 

' One of the central planks of the heavy-ion collision program is to identify and characterize the properties of 
, <~| deconfmed quark matter. Many valuable probes of the properties of the medium produced in heavy-ion collisions are 
available, that include jets, electromagnetic signals and heavy quarkonia states (QQ) In particular much work has 
\ been devoted to understand how the presence of deconfmed quarks and gluons may affect the binding of quarkonia [3] ■ 
Whether or not hadrons survive in the quark gluon plasma and how their decay widths are modified is a long- 
standing issue [|[. Since the pioneering work by Matsui and Satz [|| it is known that the screening effect due to a 
thermal medium induces the dissociation of quarkonia states. Subsequent analyses of the binding energies of mesonic 
states as a function of the temperature have been done by various authors using potential models, with an interaction 
potential extracted from lattice QCD calculations [f| |g, 0, [![. These analyses show that quarkonia dissociate at 
t N *»- ■ temperatures close to the critical temperature of QCD, T c . Analogous results have been obtained by studying the 
' correlation functions of charmonia above deconfinement Q. 

In the present paper we study the properties of heavy quarkonia states with increasing temperature. In particular we 
are interested in determining how the decay widths for the leptonic, hadronic and radiative channels arc influenced by 
the presence of the thermal medium. In order to achieve this result we first determine the radial wave function of some 
charmonia states and the corresponding binding energies. In our analysis we employ the non-rclativistic Scrhodingcr 
equation (and for comparison the relativistic Salpeter equation) with interaction potential extracted from thermal 
lattice QCD simulations. A first analysis of the decay widths of charmonia have been carried out in by studying 
the quark correlators above the deconfinment temperature. Instead we employ a wave equation and study the decay 
widths for any temperature up to the dissociation temperature of the pertinent quarkonium state. 

The remarkable result of our analysis is that some decay widths change drastically close to T c . In the transition 
region we find that the decay widths of the J/9 that depend on the value of the wave function at the origin change 
by approximately a factor 2 with respect to the corresponding values in vacuum. We also investigate the properties 
of the Xc meson and in particular of the Xcj ~ * J/9 + 7 transition. This radiative decay contributes to the total 
inclusive J/9 production by a fraction of about 0.3, as determined in proton proton and ir + ir~ data [ill ]. We find 
that this J/9 production mechanism is enhanced close to the critical temperature by approximately a factor 6 with 
respect to the value in vacuum. 

Although it has been shown in perturbation theory that the potential has a sizable imaginary part [13], we shall 
employ a real potential as extracted from lattice QCD calculations of the free-energy. Since a non pcrturbative estimate 
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of the imaginary part of the potential is not yet available, it is hard to estimate how large will be its contribution to 
the processes we are considering. 



II. THE METHOD 



Potential models have been quite successful in describing the properties of heavy quark bound states in vacuum 
(see e.g. [13]). Our key assumption is that at any temperature T, the interaction between heavy quarks can be 
approximated by an instantaneous potential, V(r,T), where r is the radial coordinate. 

We shall study bound states of heavy quarks by using the non rclativistic Schrodingcr equation 

2Mq + v{r ' T) ) ^ = mi ' (1) 

where Mq is constituent mass of the heavy quark, ipi and Ei represent the wave function and the mass of the 
corresponding QQ state, respectively. Comparing the results of this analysis with those obtained by the Salpeter 
equation for the S— wave states 



(2^/2 - V 2 + V(r, T)) fa = Miipi 



(2) 



we find excellent agreement between the outcomes of the two methods, meaning that for these mesonic states rela- 
tivistic corrections are small. 

There are some controversies about how to extract the potential V(r, T) from the color-singlet free-energy Fi(r, T) 
(see [f| 0, EH)- I n principle, in order to obtain the internal energy U\ one has to subtract the entropy contribution 
from the free-energy 



U 1 =F 1 -T 



dFi 
~dT 



(3) 



However, this subtraction procedure has been put to question, one of the reasons being that around the critical 
temperature, the potential U\ is more attractive than the potential at zero temperature. In [l4j this effect has been 
interpreted as coming from an additional interaction between the constituents in the medium. Below T c this effect 
should be due to hadrons and could be related to the string "flip-flop" interaction [l5[ ; above T c it should be related to 
the antiscreening properties of QCD. In both cases one can consider that at finite temperature the static sources Q and 
Q polarize the medium and a cloud of constituents forms around them. At short distances the main contribution to 
the potential energy comes from the QQ direct interaction, which does not change with temperature. At intermediate 
distances the polarized clouds surrounding the two heavy quarks overlap and the interaction between the constituents 
of the clouds makes the interaction potential deeper than the potential at zero temperature. This indirect interaction 
leads to a change of the interaction strength at finite temperature and is extremely strong close to T c . Some authors 
refer to this behavior of the internal energy as the "over-shooting" problem and propose to address it by defining 
V(r, T) = xUi(r, T) + (1 - x)F 1 (r, T), for < x < 1, see Ref. 0. In agreement with Ref. [3 we shall use the internal 
energy in Eq. ([3]) as the interaction potential. As we will see, with this potential we obtain dissociation temperatures 
for charmonia consistent with lattice results obtained with the Maximum Entropy Method Q . Employing the free- 
energy as interaction potential gives much smaller dissociation temperatures. Therefore, in the following we shall 
assume that x = 1; results obtained with x = are shown for comparison. 

Notice that at large distances the two heavy quarks are screened and they do not interact: in this limit the internal 
energy can be interpreted as a contribution to the constituent mass of the quark. For this reason it can be useful to 
define Ui{r,T) = U x {r,T) - E7i(oo,T). 

At short distances the static quark-antiquark free-energy is normalized in such a way that it reproduces the free- 
energy at zero temperature 



F(r,T= 0) 



a 
r 



(4) 



where a is the string tension and a is the Coulomb coupling constant. The values of these two parameters can be 
phenomcnologically fixed. Employing the numerical values a = 0.16 GeV 2 , a = 0.2 and with m c = 1.28 GeV one 



obtains M 



3.11 GeV and M v ~ 3.43 GeV. 



For non-vanishing temperatures we will consider the free-energy proposed in [If 
Dcbyc-Huckcl theory, with the parameterization of Ref. [I?) , 



employing the formalism of the 



Fi{r,T) = - 



£0/4) Vx 
2 3 / 2 r(3/4) 2 3 / 4 r(3/4) 



K\/a(x 2 + K a; 4 ) 



— (exp(— x) + x) , 
r 



(5) 
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where K1/4 is the modified Bcsscl function, x = fir, while the functions fi = (J,(T) and k = k(T) are determined by 
fitting the lattice data. Once these two functions are fixed, this parameterization of the free-energy is in excellent 
agreement with lattice data for T < 2T C . 

In Fig. [1] we compare the free-energy with the internal energy. When the free-energy is used as a potential, it 
results in an interaction that always decreases its strength with increasing temperature. Moreover the only effect that 
is present is the screening of the interaction at large distances. On the other hand U\ is more attractive close to T c than 
at T = 0, the representing curve in Fig. [1] is for T = 0.85T C . Around T c the internal energy is extremely attractive; 
further increasing the temperature it eventually becomes weaker than in vacuum, as for the case of T = 1.7 T c reported 
in Fig. m 
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FIG. 1: (color online) Full black line corresponds to the Cornell potential, Eq.Q. Full blue lines correspond to Fi and dashed 
red lines correspond to Ui, obtained with the parametrization in Eq. ([5]), for temperatures T — 0, T = 0.85 T c and T = 1.7T C . 
Notice that F\ and U\ have been normalized differently at large distances and that at vanishing temperature these two quantities 
only differ by a constant shift. 




III. MASSES AND DECAY WIDTHS 



In order to determine the decay widths it is necessary to compute the wave function and the mass of the charmonia 
states. For I = states, they are determined solving both the Scrhodinger equation and the Salpeter equation through 
the Multhopp method [l8| while for P— wave states we only solve the Schrodinger equation. For a potential with 
a singularity in r = 0, the wave function solution of the Salpeter equation is not well behaved at the origin and 
therefore we will cut-off the short distance behavior of the potential. This procedure is analogous to the one employed 
in Ref. [l9| in vacuum. In particular, we introduce a cutoff radius, r m i n , such that for r < r m i n the potential is 
constant and equal to V(r m in)- According to Ref. [2fj], we fix r m i n = In all cases we find that the typical radius 
of quarkonia states is much larger than r m i n , meaning that our results are not very sensitive to the value of the cutoff. 

Various decay processes will be considered. Annihilation processes of charmonia can be viewed as two-stage fac- 
torized process. The first process consists in cc annihilation into gauge bosons. The annihilation takes place at 
the characteristic distance r of order l/m c , so for a non-relativistic pair r — > and the annihilation amplitude is 
proportional to the wave function at the origin. Then, the produced gauge boson decays or fragments into leptons 
and hadrons. In vacuum it is assumed that the inclusive probability of the latter process is equal to one. At finite 
temperature we will assume that such a factorization persists. Therefore, the width of hadronic and leptonic decays 
will depend on temperature exclusively through the radial wave function at the origin and the mass of the meson [l3[ : 

T A oc T e+e - oc r( 3 S'i -¥ 3g) oc JJJfy ■ 

(6) 

We also consider the radiative transitions XcJ — * J + 7 from the 3 Pj levels to the 3 Si state with rate 

T B cx r XcJ _ W7 cx (2J + 1)|/ PS | 2 , (7) 



where ips is the overlap integral between the radial wave functions of the corresponding states. 
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A. A toy model for screening 



We first introduce a toy model for studying the dependence of the decay widths on the screening length. It is based 
on the assumption that at short distance the interaction potential is given by the Cornell potential and is constant 
for r > r cut : 

, n \ J -7 + or for r < r cut 

nr,r cut ) = |_^_ + (rrcut forr>rcut ■ (8) 

This procedure roughly describes what happens for non- vanishing temperatures when one uses F\ (r, T) as a potential, 
see FigfTJ We find that the \ c dissociates when r cu t — 1-04 fm, whereas the J ftp dissociates at r cu t — 0.64 fm. The 
decay widths for the processes in Eqs. ^ and (J7J) are plotted in Fig. [5] as a function of r cut . The width of the process 
XcJ J/ 1 ^ + 7 begins to decrease at r cut ~ 1.2 fm. At r cut ~ 1 fm the width for this process is zero and the Xc 
dissociates. The width for the processes in Eq. © is sensitive to smaller screening length, of less than 1 fm. These 
results are in agreement with the fact that these charmonia states have a radius of about 1 fm and that the J/<I / 
is more compact than the Xc- However, the potential in Eq. ((SJ) does not take into account the contribution from 
the "cloud-cloud" interaction [l|. As we shall see, it is not the screening of the interaction that is important in 
characterizing these decay widths, but the antiscreening of the interaction. 




FIG. 2: (color online) Charmonia decay widths as a function of r cut , obtained feeding the potential of Eq. JS} in the Schrodinger 
equation {TJ. Full blue line, decay width of Eq. (|6]); dashed red line, radiative decay width of Eq. 0. The widths have been 
normalized to the values obtained for r cut — ► 00, corresponding to the widths in vacuum. 



B. Decay widths at different temperatures 

Now we show the results for dissociation temperatures and decay widths obtained using as QQ potential the internal 
energy in Eq. ((3]) and for comparison the free-energy in Eq. ([5|) . 

We find that using the free-energy ([5]) the J/tp dissociates at 1.2 T c , while the Xc dissociate at about 0.95 T c . Using 
the internal energy as potential we find that the J/ip dissociates at approximately 1.8 T c , while the \c dissociates 
at about 1.15 T c . We have estimated the dissociation temperatures from the binding energy of the pair. Using the 
internal energy in the wave equations, one obtains higher dissociation temperatures in agreement with the fact that 
U\ is more attractive than F\. Employing the Salpeter equation we find a difference in the dissociation temperatures 
of less than 5%, in qualitative agreement with results obtained in vacuum [2l| . 

In Fig. [3]we report the values of the decay width for the process in Eq. ([6|), left panel, and for the process in Eq. {7J, 
right panel, as a function of the temperature. Employing F\ we find for both the decay processes considered here a 
behavior similar to the one reported in Fig. Therefore, one can interpret the variation of these decay widths with 
the temperature as due to the screening of the potential. On the other hand, when using U\ we find that across the 
transition region both decay widths become large. For the decays in Eq. ^ the width is about a factor 2 larger than 
in vacuum. For the radiative decay in Eq. ([7]), the ratio between the width in the thermal medium and in vacuum 
reaches a factor 6 across the transition region. This effect is not due to the screening of the potential, but is rather 
connected with the "cloud-cloud" interaction. 



5 



In order to understand whether this effect is observable one should consider the evolution of the medium in collider 
experiments. For a rough estimate, let us suppose that the \ c is produced in the early stage of the heavy-ion collision 
and then travels for about 4 fm in the thermally equilibrated medium. If the temperature is sufficiently low, so that 
the decay width is approximately the same as in vacuum, one has that less than 1% of Xc decay. On the other hand, 
for temperatures close to T c , one should observe that about 5% o_f_y c decay. Since for temperatures close to T c one can 
neglect interactions of charmonia with in-medium hadrons [H, [22j and gluons Q the radiative decay process should 
give a sizable contribution to the total decay width of the Xc- As regard the inclusive width of the J/^f, one should 
consider that at temperatures below T c , the process in Eq. ^ gives approximately the inclusive width. However, at 
larger temperatures the dominant contribution should be due to the interaction with in-medium partons 0]. 
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FIG. 3: (color online) Charmonia decay widths given in Eq. ((6]) (left panel) and in Eq. Q (right panel). We have used Fi, 
full blue lines, and Ui, dashed red lines, as potentials with the paramterization of Eq. ([5]), and normalized the widths to the 
vacuum value. 



The largest uncertainty in our calculations resides in the extraction of the potential from lattice QCD calculations 
and pertinent parametrizations to numerically evaluate the entropy. In order to test the robustness of our results one 
can employ a different parametrizations of the internal energy. Using the expression of U\(r, T) reported in Ref. Q 
we find approximately the same results reported here. In particular we find that in the crossover region the widths of 
the processes in Eqs.© and (J7J) are enhanced with respect to the vacuum values by approximately the same factors 
2 and 6, respectively. 

The analysis reported in the present paper can easily be extended to bottomia, moreover it would be of some 
interest to include D mesons. 
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